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ABSTRACT 

Wc measure the Milky Way's rotation curve over the Galactocentric range 4kpc < R < 14kpc from 
the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). 
We model the line-of-sight velocities of 3,365 stars in fourteen fields with 6 = 0° between 30° < I < 210° 
out to distances of 10 kpc using an axisymmctric kinematical model that includes a correction for the 
asymmetric drift of the warm tracer population (<tr sw 35 km s _1 ). We determine the local value of 
the circular velocity to be V c (Ro) = 218 ± 6 km s~ and find that the rotation curve is approximately 
flat with a local derivative between —3.0 km s~ kpc~ and 0.4 km s~ kpc - . We also measure the 
Sun's position and velocity in the Galactocentric rest frame, finding the distance to the Galactic 
center to be 8 kpc < Rq < 9 kpc, radial velocity Vr,q = — 10 ± 1 km s~ , and rotational velocity 
V^© = 242^3° km s -1 , in good agreement with local measurements of the Sun's radial velocity and 
with the observed proper motion of Sgr A* . We investigate various systematic uncertainties and find 
that these are limited to offsets at the percent level, ~ 2 km s _1 in V c . Marginalizing over all the 
systematics that we consider, we find that V c (Rq) < 235 km s~ 4 at > 99 % confidence. We find an 
offset between the Sun's rotational velocity and the local circular velocity of 26 ± 3 km s , which is 
larger than the locally-measured solar motion of 12 km s _ . This larger offset reconciles our value for 
V c with recent claims that V c > 240 km s _1 . Combining our results with other data, we find that the 
Milky Way's dark-halo mass within the virial radius is ^8 x 10 11 M & . 

Subject headings: Galaxy: disk — Galaxy: fundamental parameters — Galaxy: general — Galaxy: 
kinematics and dynamics — Galaxy: structure — stars: kinematics 
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1. INTRODUCTION 

The Milky Way's inner rotation curve V C (R), and in 
particular its value V c = V c (Ro) at the Sun's Galactocen- 
tric radius i?o, is crucial for our understanding of many 
Galactic and extra-galactic observations. It provides an 
important constraint on mass models for the Milky Way 
and the question of whether the Galactic disk is maximal. 
The shape of the rotation curve is an important ingredi- 
ent for realistic models of the disk's formation and evo- 
lution. The circular velocity at the Sun, which is located 
approximately 2.5 disk scale lengths from the Galactic 
center (Bovy et al. 2012b), is also important for placing 
the Milky Way in a cosmological context, for example, 
when asking whether the Milky Way follows the Tully- 
Fisher relation (e.g., Klypin, Zhao, & Somerville 2002; 
Flynn et al. 2006; Hammer et al. 2007). V c is often con- 
sidered to be an important parameter for dark-matter 
direct-detection experiments and for correcting the mo- 
tion of extra-galactic objects for the motion of the Sun, 
although we seek to dispel these notions below (§ 5.3). 

Traditionally, the local circular velocity has been ob- 
tained by measuring the Sun's motion with respect to 
an object that is assumed to be at rest with respect to 
the Galaxy. The most robust of those measurements is 
that derived from the observed proper motion of Sgr A* 
(Reid & Brunthaler 2004), even though this requires an 
estimate of i? . An alternative method that does not 
require knowledge of Ro consists of measuring the Sun's 
reflex motion using a stellar tidal stream with orbital pole 
near I = 270° (Majewski et al. 2006). Similar measure- 
ments using samples of halo stars (Sirko et al. 2004), or 
the globular cluster system (Woltjer 1975), may be con- 
taminated by the residual, presumably prograde, motion 
of these populations and, thus, only provide a lower limit. 
As discussed in detail in § 5.3, such measurements intrin- 
sically measure the Sun's rotational velocity, V^©, rather 
than the circular velocity. To arrive at V c , these mea- 
surements depend on a highly uncertain correction for 
the Sun's motion with respect to V c . The measurement 
of V c by observing the extreme linc-of-sight velocity of 
HI emission toward Galactic longitudes I ~ 90° in princi- 
ple also directly measures V^.q (Knapp et al. 1979), al- 
though the HI density drops too steeply with radius for 
V^>,0 to be directly observed, and more intricate model- 
ing is necessary to turn the HI emission toward I ~ 90° 
into a constraint on V c . 

Other measurements of V c have in the past been 
limited to measurements of the Oort constants, due 
to a lack of data at large distances from the Sun 
(e.g., Feast & Whitelock 1997). Recently, new estimates 
have been obtained from the kinematics of masers in 
the Galactic disk (Reid et al. 2009; Bovy et al. 2009a; 
McMillan & Binney 2010), and from fitting an orbit 
to the cold stellar stream GD-1 (Koposov et al. 2010). 
However, until now, no consensus has been reached as 
to whether V c lies near the IAU-rccommendcd value 
of V c = 220 km s _1 (Kerr & Lynden-Bell 1986), or 
whether it needs to be revised upward to around V c = 
250 km s" 1 (e.g., Ghez et al. 2008: Reid et al. 2009; 
Bovy et al. 2009a; Schonrich 2012). 

Most of the information relating to the shape of the 
rotation curve is based upon the observed kinematics 
of the HI emission, either through the tangent-point 



method at I < 90° and I > 270° (van de Hulst et al. 
1954; Gunn et al. 1979; Levine et al. 2008), or through 
the observed thickness of the HI layer at 90° < I < 270° 
(Merrifield 1992). Such measurements are purely geo- 
metrical and, in particular, cannot constrain the com- 
ponent of Galactic rotation that has a uniform angu- 
lar speed. This also holds for measurements of V c 
based on the line-of-sight velocities of open clusters 
(e.g., Frinchaboy 2008). The best constraints on the 
local slope of the circular-velocity curve therefore also 
come from measurements of the Oort constants (e.g., 
Feast & Whitelock 1997), although the HI observations 
are crucial to measuring the rotation curve over the full 
radial range of the disk. 

In this paper we present the first measurement 
of the Milky Way's circular velocity curve using 
kinematically-warm stellar tracers out to heliocentric 
distances of 10 kpc, from the Sloan Digital Sky Sur- 
vey Ill's Apache Point Observatory Galactic Evolution 
Experiment (SDSS-III/APOGEE; Eisenstein et al. 2011; 
S. R. Majewski, et al. 2012, in preparation). APOGEE 
is a high-resolution spectroscopic survey covering all of 
the major components of the Galaxy that — crucially — 
operates in the near-infrared, which allows stars to be 
observed to large distances in the dust-obscured regions 
of the inner Milky Way disk. While warm stellar-disk 
tracers do not on average rotate at the circular veloc- 
ity like the HI emission discussed above, their offset 
from V C (R) — the so-called asymmetric drift (Stromberg 
1946) — is a dynamical effect that can be calculated from 
their observed velocity dispersion. Because our measure- 
ment relies on a dynamical effect, we are sensitive to V c in 
the sense of the radial-force component at R, rather than 
to U^©, although the large range in I of our sample allows 
us to independently measure V^©. Additional benefits 
of using the intermediate-age to old stellar population to 
trace the dynamics of the disk are that these populations 
are much less sensitive to non-axisymmetric streaming 
motions than cold gas, and that the large asymmetric 
drift can be used to constrain the component of Galactic 
rotation that possesses a uniform angular speed. 

Using a data set of 3,365 stars at b = 0° along 14 lines 
of sight covering 30° < I < 210°, we find that the Milky 
Way's rotation curve is approximately flat over 4 kpc < 
R < 14 kpc, with V C (R ) = 218 ± 6 km s~\ A value of 
V c > 235 km s _1 is ruled out at > 99% confidence by our 
data. Our measurement of V^© = 242^° km s _1 agrees 
with the observed proper motion of Sgr A* , but the Sun's 
offset from V c is larger than the locally-measured value 
(Schonrich et al. 2010) by 14 ± 3 km s" 1 . 

The outline of this paper is as follows. In § 2, 
we describe the APOGEE data set. The methodol- 
ogy employed to model the observed kinematics of the 
kinematically-warm tracer population is presented in § 3. 
We discuss our results in § 4, including a detailed discus- 
sion of potential systematics in § 4.2. The influence of 
non-axisymmetric streaming motions on our data is as- 
sessed in § 5.1. We compare our new measurement of 
V C (R) with previous determinations in § 5.2. In § 5.3, we 
discuss the implications of our measurement of the Sun's 
offset from circular motion and, in § 5.4, we estimate 
the mass of the Milky Way implied by our data, as well 
as other recent data. We conclude in § 6. Appendix A 
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Fig. 1. — Distribution of the 3,365 mid-plane data points in 
extinction-corrected color and magnitude. The gray curves give 
approximate distances based on the peak of the distance likelihood 
in FIG. 11. 

describes how photometric distances for the stars in our 
sample are estimated. In Appendix B we discuss exten- 
sive mock-data tests used to test our methodology and 
to determine the data's sensitivity to V C (R). In what 
follows we sometimes refer to the circular velocity at the 
Sun's location V c (Rq) by the shorthand notation V c . 

2. APOGEE DATA 

The SDSS-III/ APOGEE is a near-infrared (NIR; H- 
band; 1.51 to 1.70 /im), high-resolution (R w 22,500), 
spectroscopic survey, targeting primarily red giants in 
all Galactic environments, with emphasis on the disk 
and the bulge. The APOGEE instrument (Wilson et al. 
2010, J. Wilson et al. 2012, in preparation) consists of 
a spectrograph with 300 2" fibers that reaches a signal- 
to-noise ratio of 100 per pixel (at about Nyquist sam- 
pling) at H < 12.2 in three one-hour visits during bright 
time on the 2.5-meter Sloan telescope, at the Apache 
Point Observatory in Sunspot, NM (Gunn et al. 2006). 
A detailed account of the target selection and data re- 
duction pipeline is given in G. Zasowski et al. 2012 (in 
preparation) and D. Nidever et al. 2012 (in preparation), 
respectively Here we summarize the, for our purposes, 
most important aspects of the target selection and data 
reduction. 

Our analysis is based on data from APOGEE's first 
year of regular survey operations (09-11-2011 to 05-07- 
2012). We use data from fields centered on b = 0°; 
APOGEE fields have a radius of 1.5°. By only selecting 
fields with 30° < I < 330° we avoid the bulge region. We 
resolve multiple visits to the same field by choosing the 
highest signal-to-noise ratio measurement of the Dopplcr 
shift, rather than combining the multiple epochs, be- 
cause the typical uncertainty in the line-of-sight veloc- 
ity is well below 1 km s . Only primary survey targets 



are selected, excluding targets flagged as possible cluster 
members and stars observed as part of any of various spe- 
cial programs. We exclude stars with (J — K s )o > 1.1, as 
these are problematic for the isochrone model we adopt 
to marginalize over their distances (see Appendix A). 
The resulting sample has 3,365 stars in 14 different fields 
with 30° < I < 210°. The distribution of the data 
in extinction-corrected color and magnitude is shown in 
FIG. 1 and the distribution of the stars on the sky is 
shown in FIG. 2. The properties of the sample in the 
various fields arc given in TABLE 1. 

Spectroscopic targets are selected from the 2MASS 
point-source catalog (Skrutskie et al. 2006), with the fol- 
lowing quality restrictions applied: photometric uncer- 
tainties less than 0.1 mag in J, H, and K s ; quality 
flag 'A' or 'B' in JHK S ; nearest neighbor more than 
6" away; confusion flag '0' for JHK S ; galaxy contam- 
ination flag '0'; read flag '1' or '2' for JHK S \ cxtkey 
equal to —1. The extinction corrections for targets in the 
Galactic mid-plane use mid-IR photometry from either 
WISE (Wright et al. 2010) or Spitzer-IRAC GLIMPSE- 
I (Churchwell et al. 2009); therefore, both mid-IR de- 
tections and mid-IR photometric uncertainties less than 
0.1 mag are required 24 . APOGEE's magnitude range — 
7 < H < 13.8 — is within the completeness limits for both 
the IRAC and WISE surveys. Photometry for all targets 
is extinction-corrected using the Rayleigh Jeans Color 
Excess method (RJCE; Majewski et al. 2011), which 
provides extinction values Ak with typical random un- 
certainties < 0.05 mag for individual stars using a com- 
bination of near- and mid-IR photometry. Variations in 
the adopted extinction law among different lines of sight 
can lead to differences of up to 7% (Zasowski et al. 2009). 

The color range that we select, 0.5 < (J — A' s )o < 1.1, 
includes the red clump ((J — K s )o = 0.5 to 0.7) and 
the red giant branch, for all metallicities. Dwarf con- 
tamination is most severe at (J — K s )o < 0.8, with 
most redder dwarfs being faint M and brown dwarfs. 
Most APOGEE observations consist of three or more in- 
dividual "visits" , with only three visits for stars with 
7 < H < 12.2, and six visits for 12.2 < H < 12.8. 
About 10% of our sample consists of even more visits 
for stars with 12.8 < H < 13.8. The APOGEE sampling 
is random in color, and a combination of random and 
systematic in apparent magnitude (selecting every A-th 
star in a magnitude-ordered list for each field); details of 
the selection function are unimportant for our purposes 
and it suffices to note that the spectroscopic sample is a 
representative sample of the underlying (non-extinction- 
corrected) magnitude distribution. At H < 11, dwarf 
contamination in the APOGEE fields is expected to 
be less than 10% based on population-synthesis models 
(Girardi et al. 2005) in the APOGEE fields (L. Girardi, 
2012, in preparation). Models predict that this contam- 
ination increases to about 20% at 11 < H < 12. We 
explicitly include dwarf contamination as a free parame- 
ter in the analysis below. 

Line-of-sight velocities arc determined for each individ- 
ual visit by cross-correlating against a set of « 100 syn- 

24 The first version of the APOGEE target selection did not 
consistently insist on the availability of extinction corrections, such 
that a small number of mid-plane targets do not have extinction 
estimates. We have removed 17 stars without extinction estimates 
from the sample. 
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Fig. 2. — Distribution of the data in (I, b). The smooth curve in the right panel is a uniform distribution in (I, b) over each field. 



TABLE 1 
Properties of the Sample 



Field location 

/ [degrees] 


Stars 


H < 12.2 


12.2 < H < 12.8 


12.8 < H < 13.8 


median Ak 


30° 


230 


108 


33 


89 


0.8 


60° 


111 


50 


33 


58 


0.6 


90° 


310 


178 


46 


86 


0,1 


105° 


222 


222 








0.2 


120° 


289 


157 


84 


18 


0.3 


127° 


228 


228 








0.3 


135° 


229 


229 








0.3 


142° 


227 


227 








0.6 


150° 


279 


150 


87 


42 


0.3 


165° 


225 


225 








0.1 


172° 


199 


199 








0.3 


180° 


314 


174 


90 


50 


0.2 


195° 


227 


227 








0.2 


210° 


315 


173 


91 


51 


0.2 



thetic template spectra that sparsely cover the stellar- 
parameter range 3, 500 K < T c ff < 25, 000 K in effec- 
tive temperature, T c g, —2 < [Fe/H] < 0.3 in metallic- 
ity, [Fe/H], and 2 < \ogg < 5 in surface gravity, logg 
(see D. Nidever et al. 2012, in preparation, for details 
on the exact procedure). The distribution of rms scat- 
ter in the measured line-of-sight velocity for stars with 
multiple observations peaks around 0.1km s _ . Tests of 
field-to-field variations indicate that the zero point of the 
velocity scale is stable at the 0.05 km s _1 level. A prelim- 
inary comparison between the APOGEE-measured and 
literature line-of-sight velocity of 53 stars in M3, M13, 
and M15 shows that the APOGEE zcropoint accuracy 
is w 0.26 ± 0.20 km s _1 (sec D. Nidever et al., 2012, in 
preparation for full details). In what follows we ignore 
the uncertainties on the line-of-sight velocities entirely 
(that is, we assume that they are zero), because the 
uncertainties are much smaller than any velocity differ- 
ence that we could hope to measure with our sample of 
~ 3, 000 stars having a typical dispersion of ~ 35 km s _1 . 

Repeated visits to the same star will eventually allow 
binary stars to be flagged as such based on the variability 



of their line-of-sight velocities. However, currently mul- 
tiple visits are not available for all 14 APOGEE fields 
that we employ here, and removing fields without multi- 
ple visits would significantly reduce the longitudinal cov- 
erage of our sample. Using the stars in our sample that 
have multiple visits, we find that ~4 % of stars show ve- 
locity variability at the 10 to 70 km s _1 level — the level 
at which they could be confused with the dispersion of a 
disk population. Binary contamination would spuriously 
increase the inferred velocity dispersion and therefore 
likely decrease the inferred circular velocity, due to the 
dispersion-dependence of the asymmetric-drift correction 
used below. However, as discussed below, the asymmet- 
ric drift for our sample is < 6 km s _1 at R < 10 kpc, 
such that the effect of binary contamination on the ve- 
locity dispersion is C 1 km s — . Therefore, we do not 
remove binaries from our sample, but we discuss the re- 
sults obtained when removing them in § 4.2. 

Stellar parameters and elemental abundances are de- 
termined by the APOGEE Stellar Parameters and Chem- 
ical Abundances Pipeline (ASPCAP; A. E. Garcia Perez 
et al., 2013, in preparation). After determining the 



■5 



270° 




90° 



Fig. 3. — Sum of all photometric distance distributions 
p(d\l,b, (J — K a )o, Ho, [Fe/H], DF, iso) for stars in our sample, 
viewed from the north Galactic pole, with the Sun at azimuth 
0°. These photometric-distance distributions are used in Equa- 
tion (1) to marginalize over the distance to each individual star. 
We emphasize that this is not the underlying distance distribution 
of our sample, as individual distances are poorly constrained for 
the giants in our sample. The 14 fields that we use are indicated 
by dashed white lines. 

spectral type, each star's microturbulence, T c ff, [Fc/H], 
log<?, [a/Fe], [C/Fe], and [N/Fc] arc determined through 
X 2 minimization of the difference between the observed 
and synthetic spectra derived from ATLAS9 model- 
atmosphere grids (Kurucz 1979 and more recent up- 
dates). We only use [Fc/H] in our analysis below, and 
in particular, we do not make use of logg for the pur- 
pose of selecting giants, because the measurement of 
log g has not been finalized yet within ASPCAP. Pre- 
liminary comparisons with standard stars and globular 
clusters indicate that the metallicities are currently ac- 
curate, to a precision of ^O.ldex at [Fe/H] > —1. All 
but 21 of the disk stars in our sample described below 
have higher metallicities than this, and all but 160 have 
[Fe/H] > -0.5. 

Because many stars in our sample arc on the red giant 
branch, and hence do not follow narrow magnitude-color 
relations like dwarfs and red clump stars, estimating pre- 
cise distances to these stars is difficult. In what follows, 
we avoid estimating distances to the individual stars in 
our sample by marginalizing a kincmatical model over 
the distance to the star. This marginalization is per- 
formed by integrating over the full photometric-distance 
probability distribution function (PDF) for each individ- 
ual star, obtained from models for the stellar isochrones 
and initial mass function (IMF) of the giant stars in the 
sample, as well as a prior that the stars are in an ex- 
ponential disk with a scale length of 3 kpc (Bovy et al. 
2012c). This technique is discussed in detail in Ap- 
pendix A. For dwarf stars, we assume that they are close 
enough that their distances are irrelevant for their kine- 



matics. As such, we cannot show the spatial distribution 
of the stars in our sample. As an alternative, we show in 
FIG. 3 the sum of all of the photometric distance distri- 
butions for stars in our sample (assuming that they are 
all giants). 

3. METHODOLOGY 
3.1. General Considerations 

Our approach to determining the Milky Way's rotation 
curve from the APOGEE data is to fit a kinematical, ax- 
isymmetric model to the observed line-of-sight velocities. 
In its most basic form, such a kincmatical model con- 
sists of (a) a model for the rotation curve, (b) a model 
for the distribution of peculiar velocities with respect to 
circular motion as a function of Galactoccntric radius, 
and (c) a set of parameters describing the transforma- 
tion of positions and velocities from the heliocentric to 
the Galactocentric frame. These latter parameters are 
the distance Rq from the Sun to the Galactic center, and 
the Sun's velocity with respect to the (dynamical) cen- 
ter of the Galaxy. We will ignore the small projection 
effects that arise at non-zero Galactic latitude and the 
vertical dependence of V c , since these are all at the sub- 
km s -1 level for |6| < 1.5° and distances < 10 kpc (e.g., 
Bovy & Tremaine 2012). Therefore, we arc only con- 
cerned with motions in the plane of the Galactic disk, 
and we ignore vertical motions. 

The kincmatical model provides the probability distri- 
bution of line-of-sight velocities as a function of position 
(Galactocentric radius R and azimuth <fi. or distance d 
and Galactic longitude I), after marginalizing over the 
component of the velocity tangential to the line of sight. 
This requires knowing the distance, which is only weakly 
constrained for giants (see Appendix A). Therefore, we 
additionally marginalize the kinematical model over the 
distance to each star, to obtain the probability of the 
observed line-of-sight velocity of a star, given its posi- 
tion (I, b), photometry (Jo, Hq, K s fi)i the iron abundance 
[Fe/H] (taken as representative of the overall metallicity 
of the star), the kinematical model (represented by the 
circular- velocity curve V C (R) and the distribution of pe- 
culiar velocities, DF), and the coordinate-transformation 

parameters i?o and Vq 1 = {Vr.q, V^q) (where we label 
this parameter as "gal" to emphasize that this is the full 
velocity with respect to the center of the Galaxy, not 
just the Sun's motion with respect to the local circular 
velocity) : 

p(V los \lA (J - K S ) ,H , [Fe/H], V C (R), R , V R , Q , V^, DF, iso) 
= P ( dwarf / giant ) 

dwarf /giant 

J dd P (V los \d, I, b, V C (R), R , Vr, , T^, ,DF) 

x p(d\l, b, (J — K s )o,H , [Fe/H], DF, iso, dwarf/giant). 

(1) 

The first factor within the integral is produced by 
marginalizing the kincmatical model over the veloc- 
ity component tangential to the line of sight at a 
given position. The second factor is given by the 
photometric-distance PDF, discussed in Appendix A; 
we include 'iso' in the prior information to emphasize 
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that we use an isochronc and IMF model to obtain the 
photometric-distance PDFs. We also marginalize over 
whether a star is a giant or a dwarf star, which changes 
the photometric-distance PDF; the dwarf contamination 
probability P (dwarf) is a single free parameter in our 
model. Equation (1) is the likelihood for a single star, 
and the total likelihood for the sample is calculated by 
multiplying together the individual likelihoods for all the 
3,365 data points. This total likelihood is what we opti- 
mize to fit models to the data. 

We can also calculate the distribution of line-of-sight 
velocities for individual fields, which we will use later 
to show comparisons between the best-fit model and the 
data. For the distribution of Vj os of a field centered at 
(/field, Afield) we obtain 

p{V\ os | Zfieid, Afield, V C (R), Rq, Vh i0 , V^ ,DF, iso) 

= J dldb J d(J - K s )o dHo d[Fe/H] 

p(( J - K a ) , H , [Fe/H], I, b\l Rcld , & fiold , DF, iso) 

x J ddp(V los \l,b,(J - K s ) Q ,H Q ,[Fe/R},V c (R),R , 

Vk,Q,V r 4 i 1 0,DF,iso) , 

(2) 

where p(J— K Sl H 7 [Fe/H], I, 6|/fi id, Afield, DF, iso) is taken 
directly from the field's data (that is, we sum over the 
data (l,b,{J - K s )o,H ,[Fe/E})). 

3.2. Kinematical Model 

We now describe in detail the full model that is fit 
to the data. Our fiducial model for the distribution of 
velocities in the Galactocentric rest frame is that it is 
a single biaxial Gaussian, with a mean radial velocity 
of zero because of the assumption of axisymmetry, and 
a mean rotational velocity given by the local circular 
velocity V C (R), adjusted for the local asymmetric drift 
V a (R; <tr, hn, h a ), which is a function of the velocity dis- 
persion, ctr, and the radial and radial- velocity dispersion 
scale lengths of the population (Hr and h a , respectively, 
see below). The only free parameters of the distribution 
of peculiar velocities are therefore the radial and rota- 
tional velocity dispersions, as the mixed radial-azimuthal 
moment vanishes, again because of the assumption of 
axisymmetry. The mixed moment, or equivalently the 
vertex deviation, in the solar neighborhood does not 
vanish for the warm disk population (Dehnen & Binney 
1998), but this seems likely to result from the presence 
of moving groups, and not from the smooth underlying 
distribution having a strong non-zero vertex deviation 
(Bovy et al. 2009b). 

Assuming that the distribution of velocities in the 
Galactocentric cylindrical frame is Gaussian allows us 
to analytically marginalize the velocity distribution over 
the velocity component tangential to the line of sight. 
The resulting one-dimensional distribution is itself Gaus- 
sian, with mean (V C (R) — V a {R;a,hu,h a ))sm.((f) + I) 
and variance a 2 R (l + sin 2 (0 + I) (X 2 — 1)), where X 2 is 
the ratio of the rotational and radial velocity variances 
X 2 = o-^/o-ft. The Galactocentric line-of-sight velocity 
for a star is calculated from its heliocentric line-of-sight 



velocity 



v 

Vg = Vg ISo -V R , <3 cosl + n Q R t 



!o sin / , (3) 

where we have written the Sun's velocity in terms of 
its angular velocity Hq. This angular velocity is ex- 
pected to be close to the observed proper motion of Sgr 
A* in the plane of the Galaxy (30.24 km s _1 kpc" 1 ; 
Reid & Brunthaler 2004), but we leave it as a free pa- 
rameter, in addition to Vr,©. 

The asymmetric drift in the plane of the Galaxy of a 
population of stars is given by (e.g., Binney & Tremaine 
2008) 



V c (R)V a (R) 
a 2 R (R) 



X' 



1 



d In (v*a 



d\nR 



(4) 



assuming that the radial-vertical velocity moment 
o~rz = in the Galactic plane, for symmetry reasons. 
For an exponential-disk tracer population (stellar tracer 
density i/» oc e~ R / }lR ) with an exponentially declining 
radial- velocity dispersion (ctr(R) oc e~ R / ha ), this for- 
mula becomes 



V c (R)V a (R) 1 
o-UR) 2 



X —1 + R 



1 



2 



(5) 



The main unknown in this relation is X 2 and its depen- 
dence on R. In the fit below, we assume a constant X 2 
with R for calculating the model's Gaussian line-of-sight 
velocity dispersion (see above), but to model the asym- 
metric drift more realistically, we use an X 2 {R) coming 
from an axisymmetric equilibrium model for the distri- 
bution function f(E, L) in a disk with a constant circular 
velocity with radius. For this distribution function, we 
use a Dehnen distribution function (Dehnen 1999) given 

by 



Dchi) 



i(E,L) oc 



V*{Re) 



cxp 



fiORe) [L - L C (E)] 



(6) 

where R e , L c , and Q(R e ) are the radius, angular mo- 
mentum, and angular velocity, respectively, of the circu- 
lar orbit with energy E. Using the procedure given in 
§ 3.2 of Dehnen (1999), we choose the v*(R) and <jr(R) 
functions such that they reproduce a disk with exponen- 
tial surface density and velocity-dispersion profiles to an 
accuracy of a fraction of a percent at all radii. 

The asymmetric drift for the Dehnen distribution func- 
tion (DF) for a population with <jr(Ro) = 0.2 V c (Ro) rts 
44 km s _1 is shown in FIG. 4 for hn = 3 kpc, h a = 8 kpc. 
The same is shown for a population with <7r(Ro) = 
O.lV^-Ro) w 22 km s _1 , and the difference, normalized 
for the difference in an between them, is small. For dif- 
ferent values of hn and h a we correct this function using 
Equation (5). As different values for hji and h a lead to a 
different X 2 (R), this approach is slightly incorrect, but 
tests show that the difference in V a is < 5% for reasonable 
values of Hr and h a . As the asymmetric drift is typically 
< 20 km s _1 for our sample, this leads to changes < 1 km 
s _1 , which we can ignore. Similarly, differences in A 2 for 
our sample from this model are typically < 0.2, also lead- 
ing to changes < 1 km s _1 , such that our assumption of 
a constant X 2 (R) does not strongly bias the value of the 
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Fig. 4. — Asymmetric-drift model used in this paper. The top 
panel shows the asymmetric drift V a in a Dehnen DF (Equa- 
tion [6]), expressed dimensionlessly as V c V a /o~ R . We use the 
Vr(Ro) = 0.2V C (R ) curve and show the <tr{R ) = 0.1V C (R ) 
for comparison. The middle panel shows the difference in asym- 
metric drift for different values of h R and h a , using the solid curve 
from the top panel as the fiducial model. The bottom panel shows 
the actual asymmetric-drift correction as a fraction of the circu- 
lar velocity used for the best-fit model below (solid line), and also 
shows the correction for Ro/h^ = 1 (dashed line). All curves in the 
middle and bottom panels have crfi(Ro) = 0.l4V c (Ro) (the best-fit 
value for the APOGEE sample, see below). 



inferred circular velocity. In the exponential-disk model, 
the asymmetric drift for our sample is < 25 km s _1 at 
all radii, even in the innermost disk (Freeman 1987). 

The middle panel of FIG. 4 — now for (Jr(Ro) = 
0.1414, the best-fit value below — shows the difference in 
asymmetric drift for different values of Iir and h a ; the 
main difference is at small R. This figure shows that 



the dependence of the asymmetric drift correction on h,R 
and h a at Rq is mostly < 0.02 V c — as it is unlikely that 
the scale length is as small as 2 kpc for this sample of 
stars (Bovy et al. 2012c), or that the radial-velocity dis- 
persion scale length is much shorter than 8 kpc — which 
is the same as the statistical uncertainty on the circular 
velocity that we infer below. Therefore, the systematic 
uncertainty due to the scale length of the population of 
stars is < 4 km s _1 at Ro, and likely < 2 km s _1 . We 
could directly measure the scale length of our sample, 
but as this requires a detailed understanding of the dust 
distribution in the plane of the Galaxy, and as it does 
not contribute greatly to the error budget of our mea- 
surement, we do not attempt it here. 

Our best-fit model below has a radial velocity disper- 
sion that is approximately flat over the range in R consid- 
ered here. In this case, the asymmetric-drift correction 
increases with i?, as is clear from Equation (5), and the 
difference between this dispersion profile and that consid- 
ered above is shown in the middle panel of FIG. 4. The 
bottom panel of FIG. 4 shows the actual asymmetric- 
drift correction applied in our best-fit model below. It 
is clear that this correction is small (< 6 km s" 1 at 
R < 10 kpc); hence it does not drive our inferred value 
for V c below. The bottom panel of FIG. 4 also illus- 
trates that, if the dispersion scale length were smaller 
than that in our best-fit model below, then the circular 
velocity curve would drop more steeply than what we 
infer below. 

In addition to the kinematical model discussed so far, 
we include an outlier model that consists of a Gaussian 
with a width of a = 100 km s _1 , centered on a mean 
value that is a free parameter in the model. This out- 
lier model exists primarily to deal with data-processing 
outliers, stars with a color and magnitude that places 
them beyond the disk, and close binaries with large ve- 
locity amplitudes (as we do not currently have a sufficient 
number of APOGEE epochs for each star to confidently 
remove these from the sample). In all of the fits below, 
the outlier fraction is < 1 %, thus it does not influence 
our fits. 

3.3. Dwarf Contamination 

The dwarf contribution to the likelihood in Equa- 
tion (1) could in principle be treated in the same 
way as the giant contribution, by integrating over the 
p(d\l, 6, (J — K S ) ,H , [Fe/H], DF, iso, dwarf) obtained 
from a similar isochrone and IMF model as for the giants. 
However, this is problematic, because the main sequence 
is much narrower than the giant branch, such that this 
integral is dominated by a narrow range of distances. As 
all of the dwarfs in our sample are faint main-sequence 
stars, they are expected to all be within a few 100 pc from 
the Sun, and thus be relatively unaffected by the radial 
and azimuthal gradients in velocity distribution that we 
model for the giant stars. Our approach is therefore to 
replace the integration over distance with a single eval- 
uation at zero distance from the Sun for the dwarf part 
of the likelihood. The dwarf stars are otherwise mod- 
eled using the same distribution of peculiar velocities as 
the giant stars, i.e., we assume that they are drawn from 
a similar population. The dwarf contamination in our 
best-fit model below is 7.5%, and is similar for other fits 
considered below. This is close to the value expected 
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from stcllar-population-synthcsis models (see § 2). 

3.4. Mock Data Tests 

The procedure described above makes non-trivial ap- 
proximations to an ideal axisymmetric fit to the data. 
The most important of these is that, although we include 
the asymmetric drift as part of our model, we do not in- 
clude the deviations from a Gaussian rotational- velocity 
distribution that are expected for a warm population of 
stars (with the same physical origin as the asymmetric 
drift). In addition to this assumption, we make small ap- 
proximations to the full asymmetric drift expression in 
Equation (4) . These simplifications are expected to have 
small effects on our analysis: The skew of the rotational- 
velocity distribution is only fully visible near the tangent 
point where the line-of-sight velocity is equal to the ro- 
tational velocity, but our sample covers a large range of 
distances for each line of sight, and 85% of our sample 
lies at 90° < I < 270°, where there is no tangent point. 
As such, most stars are drawn from a distribution that 
is close to Gaussian. 

To investigate these simplifications, we conduct a series 
of mock data tests. In these tests, we sample line-of-sight 
velocities from the best-fit flat-rotation-curve model be- 
low, except that we use the full Dehnen DF of Equa- 
tion (6) — uncorrected, such that it does not quite re- 
produce the assumed exponential disk — rather than the 
Gaussian approximation used in the fit (we ignore the 
best-fit X 2 as the Dehnen DF has X 2 built-in) . We then 
fit these mock data using the approximate methodology 

These mock data tests are described in detail in Ap- 
pendix B. The results from these tests show that the 
approximations made in the analysis do not produce any 
significant bias in the fitted parameters. Fits with non- 
flat rotation curves to the mock data samples indicate 
that we can reliably determine the shape of the rotation 
curve between 4 kpc < R < 14 kpc. The uncertainties 
in the fitted parameters for the mock data are approxi- 
mately the same as those for the real data (see below), 
which indicates that the uncertainties in the Galactic pa- 
rameters derived below for the real data are correct. 

3.5. Parameter Sensitivity 

We briefly discuss here how our data are sensitive to 
the parameters of the model. The mean heliocentric line- 
of-sight velocity at a position (R, (f) = (d, I) is given by 
(see Equation (3)) 

Vios = % sin (<j) + l) + Vr.q cos I - V^, sin I , (7) 

where V$ = V c — V a . We have treated the position and 
velocity of the Sun as free parameters, without imposing 
any prior on i?o, and without assuming that the solar 
rotational velocity is given by the local circular velocity 
corrected for the (previously measured) velocity of the 
Sun with respect to the local standard of rest (LSR). 
We can determine the full space motion of the Sun from 
our sample, because the longitudinal dependence of the 
correction for the solar motion of the line-of-sight velocity 
of a star at (R,(j>) = (d,l) is sinusoidal (oc V^q sin/), 
while the dependence on the rotational velocity at the 
position of the star varies differently for the d « 1 to 
10 kpc sample of APOGEE stars (oc (V c -V a ) sin(0+O)- 
Mock data tests in Appendix B show that our sample 



spans a sufficiently wide range in Galactic longitude and 
distance to disentangle these two effects, and measure the 
solar Galactocentric motion. Because we know the Sun's 
rotational frequency from assuming that the radio source 
Sgr A* at the Galactic center is at rest with respect to the 
disk, we have an independent check on our inferred value 
for the solar motion (but we do not use the measured 
proper motion of Sgr A* as a prior on the solar motion). 

We can then turn a measurement of V$ (i?) into a mea- 
surement of V C (R) by correcting for the asymmetric drift, 
V a (R), using Equation (5). This transformation requires 
a measurement, foremost, of the radial-velocity disper- 
sion, or(R). We can measure <jr almost directly using 
the I = 180° line of sight. By modeling cfr(R) as an expo- 
nential and assuming a constant tr^/crjj, we can measure 
<tr(R) from the line-of-sight-velocity dispersion at each 
I. The only subtlety here is that, in the first quadrant 
(R < Ro,l < 90°), we typically sample positions where 
Vios is close to V^, while in the second and third quad- 
rants (R > 7?o,90 < I < 270°), V\ os is mainly composed 
of Vr. We can then use the measurement of o-r(R) to 
correct for V a and derive V c . 

The fact that our data sample is composed of the 
intermediate-age to old disk population is crucial for our 
ability to measure the full shape of the rotation curve. 
Traditional tangent-point analyses of the HI emission at 
/ < 90°, or the measurement of V C (R) at 90° < I < 270° 
using the thickness of the HI layer, are insensitive to 
solid-body (uniform angular speed) contributions to the 
rotation curve, because these do not give rise to line-of- 
sight velocities for circular orbits. Formally, the HI ve- 
locities are invariant under the transformation V C (R) — > 
V C (R) + Q.R. However, this is not the case for a warm 
tracer population, as, for example, the mean rotational 
velocity transforms as 

^^V lo8 + V M (l- 1 + ^ o{R) ) rfn(* + Q. 

(8) 

where fl (R) =_V C (R)/R. Thus, the addition of QR to 
V C (R) changes V~i os by ~V a fl/^o, and rt/tt ^10% gives 
rise to a few km s _1 changes in Vi os , which we can detect 
with our data. Equation (8) clearly shows that a warmer 
population, i.e., a population with a larger asymmetric 
drift, is more sensitive to the local slope of the rotation 
curve than a colder population. 

4. RESULTS 
4.1. Basic Models 

In this section, we discuss the results from fitting the 
basic kincmatical model of a single population with ex- 
ponential i/*(R) and o~r(R) profiles to the APOGEE disk 
mid-plane data from § 2. These are the main results of 
this paper. In a subsequent section, we discuss the ef- 
fects of various systematics on these basic model fits, but 
we find in the end that they do not significantly bias the 
results for the Galactic parameters in this section. 

The results from fitting a flat rotation curve to the 
APOGEE data are given in the left column of TA- 
BLE 2. The right column of this table gives the results 
when fitting a power-law model for the rotation curve 

V C (R) = V C (R ) {R/Rof- The best-fit power-law index 



9 




30 35 40 



V Ae - V c [km s- 1 ] 



A[Fe/H] [dex] 




-. . i . . J,r. 






/ 11 v \ \ 











Multiple populations 






SFR = exp(-S/8Gyr) ■ 




Ml 





A/J [mag] 



Aj4 h [mag] 



flo [kpc] 



Fig. 5. — PDFs for the parameters of the model. All PDFs in the top two rows are for the model with a flat rotation curve, except for 
the middle two panels in the top row, which are for the power-law fit and the linear-polynomial fit. The top row has the (marginalized) 
joint PDF for V c and (from left to right) Rq, power-law index /3, linear derivative dV c /dR, and the dwarf contamination fraction P(dwarf). 



The middle row has (from left to right): the Sun's peculiar rotational velocity V, 



0,0 — ^ ne ra ti° of Ro to the radial-dispersion scale 
length, the radial-velocity dispersion, and the ratio of tangential to radial-velocity dispersion squared. The bottom row explores various 
systematics: A systematic offset in [Fc/H], in distance modulus fi, in extinction A^, and using 20 populations between 1 and 10 Gyr with 
an exponentially-declining star-formation history and an increasing velocity dispersion as in Equation (11). In the leftmost three panels in 
the bottom row, the PDF for a single systematic offset is shown as the linear density and the solid contours, while the contours of the PDF 
when allowing for a different systematic offset in inner (i < 90°) and outer (I > 90°) fields are represented as dashed lines. The figures in 
the top two rows also include the best-fit result from leave-onc-longitudinal-ficld-out fits to the data as white crosses (that is, leaving one 
of the 14 fields in TABLE 1 out of the data set). 



in this case is approximately zero, such that both mod- 
els for the rotation curve yield similar best-fit values for 
all parameters. The only significant difference between 
the two models for the rotation curve is the uncertainty 
interval for V c (Rq), which extends to lower values of V c 
in the case of the power-law fit, although the upper limit 
on V c stays approximately the same. We find that V c 
is tightly constrained in the flat rotation-curve model: 
V c = 218 ± 6 km s _1 (68% confidence). When fitting a 
power law we find V c = 218^ 9 km s _1 , with a strong 
correlation between the power-law index /3 (or, equiv- 
alcntly, the local derivative dV c /dR of the circular ve- 



locity with respect to R) and V c . Other parameters do 
not exhibit any strong correlation with V c (see discussion 
below) . 

When fitting a power-law model for the rotation curve, 
the local slope of the rotation curve is constrained to be 
-Skms-'kpc" 1 < dV c /dR < 0.4 km s" 1 kpc" 1 (68% 
confidence), with a best- fit value near the upper end 
of that range, 0.2 km s~ kpc~ . When fitting a lin- 
ear model for the circular velocity the constraints are 
similar, and we do not discuss them further; the top, 
middle panels of FIG. 5 show the joint PDF for V c 
and /3 (for the power-law fit) and dV c /dR (for the lin- 
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TABLE 2 

Results for Galactic Parameters and Tracer Properties 



Parameter 



Flat rotation curve Power-law V C (R) = V c {Ro) (R/ Ro) 



V c (Ro) [km s- 1 ] 


dV c /dR(R ) [km s _1 kpc^ 1 ] 

A [km s _1 kpc -1 ] 

B [km s" 1 kpc" 1 ] 

(B 2 — A 2 )/(2nG) [Mq pc" 3 ] 

( \ n — 1 i — In 

[km s kpc J 


218±6 
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218lt 9 

01+"'" 1 

n.2+2i 

— 2.8 
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Ro [kpc] 
Vfl, Q [km s- 1 ] 
V^, Q [km s _1 ] 
1/^,0 - V c [km s- 1 ] 

^Sgr A* I maS y r_1 l 


o 1+1.2 
8 " L -0.1 
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lu -°-0.8 
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23.915-1 

6 32+ - 07 
-° -0.70 
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D - ,3D -0.86 


a R (R ) [km s" 1 ] 
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31.4+0 1 
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Note. — Maximum-likelihood best-fit results; uncertainties are given as 68% intervals 
of the marginalized PDF for each parameter. The first block of parameters are the basic 
Galactic parameters: circular velocity V c (Ro) at the Sun's location, the derivative of the 
circular-velocity curve with R, and various other representations of these two basic parame- 
ters: the Oort parameters A and B, the rotational frequency Qq at the Sun's location, and the 
combination (B 2 — A 2 )/ (2ttG) that is relevant in the determination of the local dark matter 
density using the cylindrical Poisson equation. The second block of parameters describe the 
Sun's phase-space position in the Galaxy: distance Ro to the Galactic center, the Sun's radial 
velocity Vr q (away from the Galactic center), the Sun's Galactocentric rotational velocity 
V0 0, and the Sun's peculiar velocity with respect to local Galactic rotation V^, q — V c . We 
also give the proper motion of Sgr A* derived from the previous parameters assuming that 
Sgr A* is at rest with respect to the disk. The final block of parameters describe the tracer 
population sampled by APOGEE: radial velocity dispersion crji(R,o) at Ro, the radial scale 
length ha- of the radial velocity dispersion (given as Ro/h a ), and the ratio of the tangential 



to the radial velocity variance a^/cr^. 



ear fit). In TABLE 2, we also give the Oort constants 
A and B corresponding to the V c and dV c /di? inferred 
directly for our sample (for the case of a flat rotation 
curve, A and B are equal in magnitude and opposite 
in sign, by definition). We also present the combination 
(B 2 — A 2 )/(2irG), which is the correction term that must 
be added to the local dark-matter density to account 
for a non-flat rotation curve, when the local dark-matter 
density is inferred from local vertical kinematics using 
the cylindrical Poisson equation (e.g., Bovy & Tremaine 
2012). It is clear that we constrain this correction 
term to be smaller than a fraction of the local dark- 
matter density: -0.0025 Mope -3 < (B 2 -A 2 )/(27rG) < 
0.0004 Mq pc~ 3 , compared to a local dark-matter den- 
sity of 0.008 ± 0.003 Mgpc -3 (Bovy & Tremaine 2012). 
Thus, the correction for the non-flatness of the rotation 
curve for our best-fit parameters is approximately zero, 
and at 68% confidence it lies between — O.lGeVcm -3 
and 0.0 GeV cm -3 . We also give the rotational frequency 
fio at the Sun's location. 

Using a power-law rotation-curve model, we estimate 
the uncertainty in the rotation curve by evaluating the 
range of V C (R) at each R for 10,000 samples from the 
PDF shown in FIG. 5. Thus, every sample has a power- 
law rotation curve, but the range at each R does not have 
to follow a power law itself. This exercise shows whether, 
in the power-law model, the value of V C (R) is well- 



constrained or not for each R; The result is shown in the 
left panel of FIG. 6. The same is shown in the right panel 
of that figure for a fit assuming a cubic-polynomial model 
for the rotation curve. However, in the latter model we 
have imposed a prior that Rq < 9 kpc, as a large number 
of samples from the PDF have Rq > 9 kpc, which we as- 
sume to be implausible from prior data (e.g., Ghez et al. 
2008; Gillessen et al. 2009). We emphasize that we do 
not include such a prior in any other part of the analy- 
sis presented here. We see that the rotation curve is best 
constrained at small radii (4 kpc < R < 8 kpc), and that 
it is largely consistent with being close to flat over the 
full range of Galactocentric radii considered here. 

We also determine the full planar phase-space position 
of the Sun in the Galactocentric reference frame. The 
current APOGEE data arc relatively insensitive to the 
value of Rq. The 68% interval for Rq for both the flat- 
rotation-curve and power-law fits is about 8 kpc < Rq < 
9 kpc, which is similar to that found for the mock data 
in Appendix B. The best- fit value is at the lower end of 
this range: 8.1 kpc when assuming a flat rotation curve, 
8 kpc for the power-law rotation-curve fit. The fully 
marginalized PDF for R in the top, left panel of FIG. 5 
is almost entirely flat between 8.0 kpc < Rq < 9.3 kpc. 

The Galactocentric motion of the Sun is approx- 
imately Vr,© = —10.5 ± 1.0 km s^ 1 , and = 
242I3 km s~ (in the case of a flat rotation curve) or 
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Fig. 6. — The Milky Way's rotation curve in the range 4 < R. < 14 
kpc as inferred from our data for different forms of the shape of the 
rotation curve. At each R, the range in V C (R) from 10,000 samples 
of the PDF — assuming a power-law or cubic-polynomial model for 
the shape of the rotation curve — is determined and the 68%, 95%, 
and 99% intervals are shown in varying shades of gray. For the 
cubic-polynomial model, we impose a prior that Rq < 9 kpc (see 
text). For comparison, the squares show the rotation curve of M31 
from the compilation of Carignan et al. (2006). 



V^© = 242^17 km s 1 (for a power-law fit to the rota- 
tion curve). In the latter fit, there is a strong correlation 
between V c and V^©; the difference between the two is 

23. l±li km s -1 . 
V c in FIG.' 5 is well- 
-l 



much better constrained: V c 



<t>,Q 



V, 



The marginalized PDF for V^© 

described by V^© — V c = 26 ± 3 km s _i . We discuss the 
consequences of this solar motion in detail in § 5.3, but 
we note here that the estimate of the angular motion of 
the Galactic center that we obtain from combining our 
fits for V^.q and i?o is consistent with the proper motion 
of Sgr A* as measured by Reid & Brunthaler (2004): our 
estimate is /xs gr a* = 6.3lg'ymas yr _1 , compared to the 
direct measurement of 6.379 ± 0.024 mas yr _1 . We dis- 
cuss the apparent discrepancy between our agreement 
with the proper motion of Sgr A* and our low value for 
V c in § 5.3. 

The final block of parameters in TABLE 2 describe 
the tracer population. The velocity dispersion that we 
infer for the tracer stars is close to that expected for 
an old disk population: (7r(Ro) ~ S2.0^ 5 km s _1 for 
the flat-rotation-curve and power-law fits. The ratio of 
the tangcntial-to-radial velocity dispersions squared is 
0.69 < X 2 < 1.0, with the best-fit value at the lower end 



of this range. This value is higher than expected from 
the epicycle approximation for a flat or falling rotation 
curve, which is X 2 < 0.5. However, this expectation 
holds only for a cold disk, and corrections due to the 
temperature of the old disk population always increase 
X 2 near R (Kuijken & Tremaine 1991): the Dehnen 
disk distribution functions of Equation (6) have X 2 that 
varies spatially, and reaches approximately 0.65 near Rq 
(Dehnen 1999). The best-fit value for Ro/h a is approx- 
imately zero, with non-zero positive values ruled out by 
the data: the 68% interval is -0.24 < R /h a < 0.03. 
Thus, the radial- velocity dispersion does not drop expo- 
nentially with radius with a scale length between 2 Rq/3 
and Rq', such a drop would be expected from previ- 
ous measurements of the radial dispersion as a func- 
tion of R (Lewis & Freeman 1989), or from the observed 
exponential decline of the vertical velocity dispersion 
(Bovy et al. 2012b) combined with the assumption of 
constant a z /aR. Wc have attempted fits with two popu- 
lations of stars with different radial scale lengths (3 kpc 
and 5 or 6 kpc) and radial-velocity dispersions, but the 
same radial-dispersion scale length. The best-fit Ro/h a 
remains zero, such that it does not seem that we are see- 
ing a mix of multiple populations that conspire to form 
a flat <7r profile. 

Even with the best-fit flat radial-dispersion profile, the 
disk is stable over most of the range in R considered 
here. The Toomre Q parameter — Q = an k/(3.36GE) 
(Toomre 1964) — for a flat rotation curve is 



Q = 1.72 



°R 



32 km s 



R 



8 kpc 



220 km s 
S 

50 M© pc" 2 



(9) 



This expression has Q > 1 down to 4.9 kpc and Q = 0.91 
at R = 4 kpc for a constant (Jr(R) and a surface den- 



sity S oc e 



-R/(3 kpc) 



Although the disk is marginally 



unstable in our best-fit model, this conclusion depends 
strongly on the assumed radial scale length: for Iir = 
3.25 kpc, Q > 1 everywhere at R > 4 kpc. The flat- 
ness of the inferred or profile also depends on the as- 
sumed constancy of X 2 . Actual equilibrium axisymmet- 
ric disks, such as those having a Dehnen distribution 
function (Equation (6)), have a radially-dependent X 2 , 
with X 2 at R = 4 kpc typically smaller than at R = 8 to 
16 kpc (Dehnen 1999, Figure 4). At R = 4 kpc, which for 
the present data sample is only reached for the I = 30° 
line of sight, the line-of-sight velocity is entirely com- 
posed of the tangential- velocity component, such that 
any decrease in X 2 leads to an increase in an to sustain 
the same <j$. Therefore, the true an(R) profile presum- 
ably is falling with R, and the entire disk at R > 4 kpc 
should be stable in our model. 

Full PDFs for all of the parameters of the basic models 
discussed in this section are given in FIG. 5. It is clear 
that with the exception of V c and the derivative of the 
rotation curve — j3 in the power-law model and dV c /di? 
in the linear model — there are no strong degeneracies 
among the parameters. Also included in this figure are 
the results from fitting all but one of the 14 APOGEE 
field for each field: these leave-one-out results show that 
no single field drives the analysis for any of the paramc- 
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ters. 

We compare the best-fit flat-rotation curve model with 
the data in FIG.s 7 and 8. FIG. 7 shows the raw he- 
liocentric velocities versus Galactic longitude for all of 
the stars in the sample, as well as the mean for each 
field, and the prediction for the mean from the best-fit 
model. This "model mean" is computed by construct- 
ing the predicted line-of-sight velocity distribution us- 
ing Equation (2) and then calculating its mean. The 
full predicted line-of-sight velocity distribution for each 
field is shown in FIG. 8 as the smooth curve, which is 
to be compared with the histogrammed data points for 
each field. It is clear from FIG. 8 that the fits capture 
the overall smooth features of the distribution of line-of- 
sight velocities in each field, including any asymmetry in 
the distribution. However, it is also obvious that some 
of the model predictions are shifted from the observed 
histograms. This effect also appears in the comparison 
between the data mean and model mean for each field in 
the lower panel of FIG. 7. There appears to be a distinct 
pattern in these residuals. Thus, the best-fit axisymmet- 
ric model fails to account for w 10 km s _1 smooth devi- 
ations in the velocity field. These deviations are likely 
due to non-axisymmetric streaming motions, which are 
expected at this level of accuracy; We discuss this fur- 
ther in § 5.3. With the exception of the I = 60° field, 
the data seem to follow the smooth velocity distribution 
quite well, albeit somewhat shifted on average. 

4.2. Systematica 

The basic models presented in the previous section are 
potentially subject to systematic uncertainties, which we 
discuss in this section. These uncertainties are mainly 
related to the distances of the stars in our sample, as 
these are obtained from somewhat uncertain models of 
the color-magnitude distribution of giant stars in the 
near-infrared J, H, and K s bands. We discuss here the 
influence of our assumed value for the scale length of the 
stars, of potential systematic offsets in the preliminary 
APOGEE metallicitics used to inform the photometric 
distances, of systematic distance uncertainties, and of 
the impact of mis-estimated interstellar extinction and 
absorption. We also discuss the effect of binary con- 
tamination, of using a different set of stellar isochrones 
for distance estimation, and of modeling the data using 
many underlying populations of stars, with velocity dis- 
persions smoothly varying with age. 

In our fits, we have simply assumed a value of 3 kpc for 
the radial scale length, hn, of the sample of stars used in 
the analysis. This assumption is based on the fact that, 
locally, the dominant solar-mctallicity disk population 
has this scale length (Bovy et al. 2012c). The parameter 
hji enters our analysis in two ways, the first being in the 
asymmetric-drift correction (Equation (5)), and the sec- 
ond being in the prior on the distance distribution (Equa- 
tion (Al)). The former of these is the most important 
for determining V C (R). Changing ha to 2 kpc — which is 
highly unlikely, given that such a short scale length is 
not observed for any of the metal-rich, "thin disk" pop- 
ulations in the solar neighborhood (Bovy et al. 2012c) — 
increases the best-fit value to V c = 223 km s - for a 
model with a flat rotation curve, and to V c = 220 km s _1 
for a power-law rotation curve. Increasing hn to 4 kpc 



decreases the best-fit value to V c {Rq) = 217 km s~ . 
Similar changes happen for the mock data sets described 
in Appendix B. Therefore, we conclude that the assumed 
radial scale length has a negligible influence on the in- 
ferred V c . 

To determine the influence of systematic offsets in the 
metallicities used in the photometric distances, we have 
performed fits allowing for a single metallicity offset, 
A[Fe/H], for all of the stars in the sample, for a model 
with a flat rotation curve. We find that the best-fit 
A[Fe/H] = -0.15 ± 0.02 dcx, but the best-fit value of 
V c is unchanged, as is its uncertainty (see FIG. 5 for 
the A[Fe/H]-V c PDF). Although an offset of -0.15 dex 
might seem large, this offset should not be thought of as 
a measurement of [Fe/H] for the stars in our sample, 
because the photometric-distance PDFs, such as that 
in FIG. 11, arc relatively insensitive to large changes 
in [Fe/H], especially at the high-mctallicity end. We 
have also performed a fit with a different A [Fe/H] for 
inner (I < 90°) and outer-disk fields. Such a fit finds 
A [Fe/H] = -0.12 ± 0.03 dex for the outer-disk fields, 
and A [Fe/H] = 0.70 ± 0.06 dex for the inner-disk fields, 
with a best-fit V c = 220 km s" 1 ; see FIG. 5 for the full 
PDF of the A [Fe/H] and V c . As discussed previously, 
the photometric-distance PDFs arc largely insensitive to 
changes in [Fe/H] at high [Fc/H]; the Padova isochrones 
also have an upper limit of [Fc/H] max = 0.45 dcx, such 
that increasing a star's [Fe/H] beyond this has no effect, 
because we then use [Fc/H] max . Therefore, even a change 
of A[Fc/H] = 0.7 dex has a negligible influence on the in- 
ner, high-metallicity, disk fields. In none of these fits is 
our inferred range for Rq affected. 

Similarly, we have performed fits allowing for a sin- 
gle distance-modulus offset, A/j,, to be applied to all 
of the stars in our sample, or for a separate offset for 
inner- and outer-disk fields. In the former case, we find 
A/i = -0.27 ± 0.06, corresponding to a 12 ± 3% dis- 
tance offset, without any influence on the inferred value 
of V c - In the inner/outer offset fit, we find 23% smaller 
distances at I < 90°, and 20% larger distances in the 
outer disk, again with a negligible influence on the best- 
fit V c (V c = 221 km s^ 1 ), but with a second minimum 
around V c = 190 km s^ , and a wider V c PDF that is, 
however, at V c < 235 km s" 1 at 99% confidence. A more 
physical effect is to allow for an offset in the extinction, 
A Ah, again using a single offset for the whole sample, 
or splitting the sample into inner- and outer-disk fields. 
A single extinction offset is well-constrained to be small: 
AA H = -0.01 ± 0.01. Similarly, the best-fit extinction 
offset in the outer disk is AAh = ± 0.02, while in the 
inner disk we find some evidence for mis-estimated ex- 
tinction, with AAh = — 0.2±0.03. In both of these cases, 
the best-fit V c increases by 1 km s _1 or less, and the un- 
certainties remain approximately the same (see FIG. 5). 

As discussed in § 2, we do not limit the sample to stars 
for which multiple velocity epochs are available that in- 
dicate that they are not part of a multiple system. Such 
a cut significantly reduces the longitudinal coverage of 
our sample. When we do apply this cut, we find that the 
best fits for V c and all other parameters are unchanged 
from the results for our basic flat-rotation-curve model in 
§4.1, albeit with larger uncertainties. We find in partic- 
ular that V c = 217 ± 7 km s~ , but the uncertainties in 
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Fig. 7. — Comparison between our best-fit flat-rotation-curve model and the data. The top panel shows the individual heliocentric 
velocities of stars in our sample as a function of Galactic longitude, with the mean for each APOGEE field indicated by a gray dot. The 
prediction for this mean from the best-fit model is shown as the gray cross. The bottom panel shows the difference between the data mean 
Klata and the model mean V mo d e i. Note that the residual for the I = 60° field in the bottom panel is missing as it is larger than 15 km s _1 . 



the position and velocity of the Sun are much increased. 

In Appendix A. we describe how we employ Padova 
isochrones to estimate photometric distances to the stars 
in our sample. When instead we use isochrones from 
the BaSTI library (Pietrinferni et al. 2004), using the fil- 
ter transformation from Carpenter (2001) to transform 
K s to the if-band filter used by BaSTI, we find that 
the results of our basic models (both with a flat and 
a power-law rotation curve) are essentially unchanged 
(with changes in V c of about 1 km s~ ). 

To determine whether our results change when we fit 
the data as a mix of multiple stellar populations, with 
velocity dispersions varying from that of a young popu- 
lation to that of the oldest disk population, we have per- 
formed a fit where the likelihood for Vi os in Equation (1) 
is generalized to 

P (V los \l, 6, (J - K s )o, H , [Fe/H], V c (R),Ro,Vr, @ , V^ J0 ,iso) 
= Y, p ( l ) P(Vios\l, b, ( J - K s ) , H , [Fc/H], V C (R),R , 

i 

V R , Q ,V^ e ,DFi,iso), 

(10) 

where DF^ is the single- Gaussian population model used 
before, with dispersion at age r. 



n +0.1 Gyr 
10.1 Gyr 



0.38 



(11) 



and using an exponentially-declining star-formation his- 
tory, such that P(i) oc e Tl/(8Gyr) . This model follows 



the fit of Aumer & Binney (2009) in the solar neighbor- 
hood. We use a mix of 20 equally age-spaced populations 
between 1 and 10 Gyr, and find V c = 216 ± 7 km s _1 , 
showing that this alternative model does not change our 
results; the full V c -R PDF is shown in FIG. 5. This 
model has crj^o = 38 ±2 km s -1 , as expected for the old, 
"thin disk" population. 

5. DISCUSSION 
5.1. Effect of Non-Axisymmetric Features 

In the fits in § 4, we have assumed axisymmet- 
ric models for the velocity field in the Galactic disk. 
However, the Milky Way is known to have some 
non-axisymmetric structures, such as the central bar 
(Blitz & Spergel 1991; Binney et al. 1991), spiral arms 
(e.g., Drimmel & Spergel 2001), and a potentially triax- 
ial halo (Law et al. 2009), that may influence the mea- 
surement performed in this paper. Because we use a 
warm disk population as dynamical tracers, we can ex- 
pect the influence of any non-axisymmctry to be smaller 
than for the colder gas tracers used in other studies (e.g., 
Levine et al. 2008), because warm tracers respond less 
strongly to non-axisymmetric perturbations to the Milky 
Way potential than cold tracers (Lin et al. 1969). 

We have calculated the velocity field for a warm 
tracer population in various non-axisymmetric models 
(see J. Bovy, 2013, in preparation, for full details). We 
follow the procedure of Dehnen (2000) and Bovy (2010) 
to calculate the velocity distribution at a given position 
for a given non-axisymmetric model for the disk, by per- 
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Fig. 8. — Comparison between the full line-of-sight heliocentric 
model. Fields are ordered by increasing longitude. The agreement 
offsets in the mean of the distributions, as is also clear from FIG. 7. 



velocity distributions of the data and the best-fit flat-rotation-curve 
between the gray model curve and the data is good, except for some 



forming backward orbit integrations until before the on- 
set of the non-axisymmetric perturbation in an axisym- 
metric disk. At that time, the distribution function of 
the pre-existing axisymmetric disk is evaluated. By the 
conservation of phase-space density, this probability is 
equal to the current probability of the phase-space start- 
ing point of the orbit integration. We assume a fiat rota- 
tion curve for the axisymmetric part of the potential, and 
we model the initial, axisymmetric distribution function 



as a Dchncn distribution function (Equation (6)) with a 
velocity dispersion of 0.2 V c , a radial scale length of Rq/3 
and a dispersion scale length of Rq. These parameters are 
close to the best-fit parameters for the APOGEE tracer 
population. 

We calculate the mean velocity field for the bar model 
of Dehnen (2000), and use it to construct the mean, 
non-axisymmetric Vi os field. This field is then applied 
as part of the model, and the result for the mean line- 
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Fig. 9. — Influence of non-axisymmetry on the predicted mean of the Vi os distributions. The flat-rotation curve model of TABLE 2 and 
FIG.s 7 and 8 is displayed in the uppermost panel. Non-axisymmetric models shown are: the model for the Milky Way bar of Dehnen (2000); 
a two-armed logarithmic spiral with a fractional amplitude of 1% and a pitch angle of —15°; elliptical disk models of Kuijken & Tremainc 
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of bar and spiral structure perturbations is small, an elliptical disk model can reduce the size of the residuals. This figure only shows the 
influence of non-axisymmctry on the mean of the line-of-sight velocities in each field; much more information is of course contained in the 
full distribution of linc-of-sight velocities. 



of-sight velocity in each field is shown in FIG. 9 — there 
is of course much more discretionary power in the full 
distribution functions. In the same manner, we calcu- 
late the mean V| os field for a two-armed logarithmic spi- 
ral with a fractional potential-amplitude of 1% of the 
background, axisymmetric potential, a pitch angle of 
— 15°, an angular frequency of 0.65 fio (placing the Sun 
near the 4:1 inner Lindblad resonance), and an angle 
between the Sun-Galactic-center line and the line con- 
necting the peak of the spiral pattern at the solar radius 
of 20°. We do the same for models with a flat ellipti- 
cal (stationary, m = 2) distortion to the potential (see 
Kuijken & Tremaine 1994), with a fractional amplitude 
of 5% and position angles of 0° and —45° with respect 
the the Sun-Galactic center line. All of these perturba- 
tions were adiabatically grown. We see in FIG. 9 that 
the influence of the bar and spiral structure is negligi- 
ble in the mean Vi os velocity field. An elliptical, m = 2, 
perturbation could distort the mean velocity field in a 
way that would affect our analysis. Naively, the model 
in the bottom panel of FIG. 9_ significantly reduces the 
residuals between the best-fit V\ os and the data Mos- 
It is clear that the data used in this paper could be 



used to constrain non-axisymmctric perturbations to the 
axisymmctric potential assumed here. We defer a full 
treatment of this to a subsequent paper. 

5.2. Comparison with Other Determinations of the 
Local Circular Velocity 

There have been many previous determinations of the 
local circular velocity, and we discuss here how our new 
measurement compares to these. Determinations of the 
local circular velocity can be roughly divided into two 
groups: measurements of the Sun's velocity with respect 
to an object or population assumed to be at rest with re- 
spect to the Galactic center (e.g., Sgr A* or a population 
of halo objects), or direct measurements of the local ra- 
dial force (e.g., by determining the Oort constants or the 
orbit of a stream of stars) . The former directly measure 
the Sun's Galactocentric velocity, V^q, but must assume 
a value for the Sun's motion with respect to the circular 
orbit at Ro to arrive at V c (Ro). The most direct mea- 
surement of this kind is the combination of the precisely 
measured proper motion of Sgr A* (Reid & Brunthaler 
2004) with the distance to the Galactic center deter- 
mined from the Keplerian orbits of S stars in the in- 
nermost parsec (Ghez et al. 2008; Gillessen et al. 2009). 
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These measurements give i?o ~ 8 kpc, which combined 
with the proper motion of Sgr A* of 30.24 km s _1 kpc" 1 . 
yields a total solar velocity of 242 km s . As noted be- 
fore, this result is entirely consistent with our inferred 
value of the angular motion of the Galactic center, and 
with our Vfao = 242^° km s _1 . The discrepancy be- 
tween the value of V c = 229 km s" 1 in Ghez et al. (2008) 
(or higher measurements also based on the proper mo- 
tion of Sgr A*) and our V c = 218 km s" 1 is therefore 
entirely due to our different value for the solar velocity 
with respect to V c , and intrinsically these measurements 
are consistent. Similarly, the recent measurements of 
V^ Q = 276 ±23 km s" 1 and F^ >0 = 244 ± 14 km s" 1 us- 
ing the kinematics of the Sgr stream (Carlin et al. 2012) 
are consistent with our measurement of the solar velocity. 

Measurements based on the Oort constants (e.g., 
Feast & Whitelock 1997), or the dynamics of a cold stel- 
lar stream (Koposov et al. 2010), also indicate that V c ~ 

220 km s~ . The measurement of the Oort constants 
from Hipparcos proper motions by Feast & Whitelock 
(1997) directly measure the angular frequency of the local 
circular orbit, independent (in principle) from the solar 
motion: O = V c /R = 27.2 ± 0.9 km s _1 kpc _1 . This 
value is consistent with our determination (see TABLE 2) 
Q a = 27.0±3 j km s" 1 kpc" 1 . The determination of V c = 

221 ± 18 km s" 1 of Koposov et al. (2010) from fitting an 
orbit to the cold GD-1 stellar stream at R ~ 10 kpc is 
clearly consistent with our measurement, although note 
that they assumed the V^© — V c = 5.25 km s" 1 value for 
the solar motion from Dehnen & Binney (1998). 

The recent measurement of V c = 254 ± 16 km s~ , 
from the kinematics of mascrs in the Galactic disk 
(Reid et al. 2009), was shown to have an overly opti- 
mistic precision by Bovy et al. (2009a), who concluded 
from a more general model for the distribution func- 
tion of the masers that these data only imply V c = 
246 ±30 km s" 1 . The measurement of V c from the maser 
kinematics also assumes that V^.q is given by V c plus the 

locally-determined solar motion of ~ 10 km s~ . Apart 
from indicating a high value for V c , albeit with a large 
uncertainty, the masers were also found to be lagging 
with respect to circular motion by about 15 km s~ , a 
large offset for a young and relatively cold tracer pop- 
ulation. This offset cannot be explained by the asym- 
metric drift and it is physically implausible since it re- 
quires the masers to have more eccentric orbits than most 
of the young stars (McMillan & Binney 2010). Requir- 
ing the masers to be on circular orbits assuming a flat 
rotation curve leads to a solar motion of V^q — V c = 
18.6±2.4 km s" 1 and V c = 232±24 (McMillan & Binney 
2010), in good agreement with our measurements. Our 
measurement, however, relies on the theoretically well- 
motivated asymmetric-drift correction that is a direct 
consequence of the collisionless Boltzmann equation. 

Therefore, we conclude that our measurement of V c is 
consistent with previously reported values. Compared 
to other determinations, our measurement of V c is one of 
the most highly precise, and unlike many other measure- 
ments, it does not require assumptions about the relation 
between the solar motion and V c . 

Our data strongly rule out that the Milky Way's 
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Fig. 10. — Sun's peculiar velocity with respect to circular motion. 
This figure shows the posterior probability distribution function of 
the Sun's velocity with respect to the circular orbit at the solar 
radius. The measurement of the Sun's motion with respect to the 
orbit of a zero velocity dispersion population from Schonrich et al. 
(2010; SBD10) is indicated by dashed lines. The fact that our mea- 
surement of the Sun's rotational peculiar velocity does not agree 
with the SBD10 value may indicate that the closed orbit at the 
solar radius is not circular, due to non-axisymmetry at the level 
of 10 km s _1 . The inset shows this non-circular closed orbit (spa- 
tially to scale) if the non-axisymmetry is due to ellipticity of the 
disk (e.g., Kuijken & Tremaine 1994); the dashed gray curve is the 
circular orbit. Alternatively, the locally-determined solar motion 
may be off by ~ 10 km s . 



circular velocity at the solar radius is > 235 km s 
Marginalizing over all of the systematics discussed in 
§ 4.2, V c < 235 km s -1 at > 99% confidence. Fixing 
V c to 250 km s _1 with Rq = 8 kpc or 8.4 kpc leads to 
a best- fit with much larger % 2 : A\ 2 = 47 (~ 5er) and 
Ax 2 = 34 (~3.7er), respectively. 

Our measurement that the Milky Way's rotation 
curve is essentially flat over 4 kpc < R < 14 kpc 
agrees with measurements based on the kinematics of 
HI emission (e.g., Gunn et al. 1979; Merrifield 1992) 
and with the local measurement of the Oort constants 
(Feast & Whitelock 1997). 

5.3. Implications for the Motion of the LSR 

We found in § 4, in both the flat- and power-law- 
rotation-curve models, that the Sun's velocity with re- 
spect to the center of the Galaxy — a distinct param- 
eter from V c in our fit — is larger than V c (Rq) by ~ 
24 km s . Defining the Rotational Standard of Rest 
(RSR; Shuter 1982) as the circular velocity in the ax- 
isymmetric approximation to the full potential, this solar 
motion with respect to the RSR is much larger than that 
measured by applying Strombcrg's asymmetric-drift re- 
lation to local samples of stars. The asymmetric-drift 
relation is used to estimate the velocity of the zero- 
dispersion orbit by extrapolating the relation given in 
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Equation (5) from warmer samples of stars to or = 0; 
it is this zero-dispersion orbit that is typically what 
is meant by the LSR (Fich & Tremaine 1991). Such 
analyses of the Hipparcos data yield a solar motion 
with respect to the LSR that is somewhere between 
5 and 13 km s _1 (Dchnen & Binney 1998; Hogg et al. 
2005; Schonrich et al. 2010); 13 km s~ is also the offset 
of the average rotational velocity of nearby stars (e.g., 
Allende Prieto et al. 2004). The discrepancy between 
our determination of the solar velocity with respect to 
the RSR, and that with respect to the LSR, is clearly 
shown in FIG. 10, where the PDF for V^ a - V c , Vr.q is 
shown for the flat-rotation curve model, together with 
the locally-measured motion with respect to the LSR 
of Schonrich et al. (2010). The analogous hgure for the 
power-law rotation-curve fit has an essentially identical 
appearance. It is clear that the locally-measured solar 
motion V^,q is not equal to our globally-measured solar 
motion to very high significance. Or, equivalently, the 
motion of the RSR is not the same as that of the LSR: 
the LSR seems to rotate ~ 12 km s _1 faster than the 
RSR. 

This discrepancy may result from a breakdown in 
the assumptions of the locally-measured solar motion 
using Strdmberg's asymmetric-drift relation. As con- 
vincingly shown by Schonrich et al. (2010), the neglect 
of the radial metallicity gradient in the analysis of 
Dehnen & Binney (1998) leads to a correction of <~ 
7kms~ 1 , but further improvements in the chemical- 
evolution model of Schonrich et al. (2010) might lead 
to further corrections. The local velocity distribution 
has also long been known to contain various streams 
or "moving groups" (e.g., Kapteyn 1905; Schwarzschild 
1907; Eggen 1986; Dehnen 1998; Bovy et al. 2009b) that 
are likely of a dynamical origin (Bensby et al. 2007; 
Famaey et al. 2008; Bovy & Hogg 2010; Sellwood 2010) 
and contain a significant fraction of the nearby stars 
(Bovy et al. 2009b). These moving groups make an ac- 
curate determination of the solar motion challenging. As 
such, the locally-determined correction for the Sun's mo- 
tion with respect to the LSR may well be incorrect by 
~ 5 km s~ or more. A solar motion of ~ 24 km s - 
rather than ~ 12 km s _ is not much more unlikely given 
the distribution function for a ~ 5 Gyr population: using 
a Dehnen distribution function (Equation (6)) with a ra- 
dial dispersion of 30 km s _1 (Equation (11)) and assum- 
ing our best-fit model gives P(V^© — V c > 12 km s _1 ) w 
0.2 and P(V^ G - V c > 24 km s _1 ) « 0.1. 

However, if Strdmberg's relation leads to a solar mo- 
tion with respect to the LSR that is different from the 
globally-determined solar motion with respect to the 
RSR, the most straight-forward interpretation is that 
LSR's orbit is not circular, but deviates from circu- 
larity by about 10 km s~ x because of large-scale, non- 
axisymmetric streaming motions. The deviation cannot 
be much larger than this value, as the radial compo- 
nent of the solar motion as determined from local sam- 
ples agrees with our global measurement. As our mea- 
surement of the circular velocity is mainly at azimuths 
(j) < 45°, this 10 km s _1 deviation must happen over this 
angular scale. 

There is mounting evidence for the existence of stream- 



ing motions in the Galactic disk of this magnitude. There 
has been a long-standing discrepancy of ^7 km s~ 1 be- 
tween the rotation curves determined from first (I < 90°) 
and fourth (I > 270°) quadrant tangent-point measure- 
ments (e.g., Kerr 1962; Gunn et al. 1979; Levine et al. 
2008), which could indicate streaming motions. Simi- 
larly, the analysis of the extreme line-of-sight velocity 
of HI emission toward I ~ 90° leads to differences in V c 
of approximately 30 km s _1 , albeit with large uncertain- 
ties (Knapp et al. 1979; Jackson & Kerr 1981; Jackson 
1985). The fact that the maximum linc-of-sight velocity 
in both CO and 21-cm data reaches zero at I ~ 70°, rather 
than at I = 90°, has been used to argue that the LSR 
moves ahead of V c with a speed of ~7.5 km s _1 (Shuter 
1982; Clemens 1985), similar to the offset we find. More 
recently, an analysis of line-of-sight velocities from the 
RAVE has found a gradient, dVfj/di?, in the mean radial 
velocity, Vr, of ~ 3 km s~ kpc - , and large streaming 
motions in both Vr and the mean tangential velocity, V^ , 
within a few kpc from the Sun (Siebert et al. 2011). 

As a simple exploration of this possibility, we have 
computed the closed orbit at Ro in a model for the Milky 
Way disk where this streaming motion is due to elliptic- 
ity of the disk, following Kuijken & Tremaine (1994). We 
use a flat rotation curve, with a cos(2 <f>) perturbation of 
constant ellipticity having an amplitude of 14 km s _1 . 
Such a perturbation could arise if, for example, the halo 
is triaxial (e.g., Law et al. 2009). The closed orbit is 
shown to scale in the inset in FIG. 10; it has an eccen- 
tricity of 0.06. 

Regardless of the origin of the discrepancy between our 
globally-measured solar motion and the locally-measured 
value, our larger preferred value of V^© — V c means 
that the Sun is likely closer to the pericenter of its or- 
bit around the Galactic center than previously believed. 
Revising V^,© — V c upward by 12 km s _1 increases the 
eccentricity of the Sun's orbit to ~ 0.1 from ~ 0.06 and in- 
creases its mean Galactocentric radius to ^8.9 kpc from 
~ 8.5 kpc. Such a large mean radius increases the ten- 
sion between the Sun's high metallicity compared to local 
stars (Wielen et al. 1996), although the peak of the local 
metallicity distribution may be closer to solar metallicity 
than previously believed (e.g., Casagrande et al. 2011). 

Finally, we address what our results imply for dark- 
matter direct-detection experiments, and for correcting 
the motion of Galactic and extra-galactic objects for the 
motion of the Sun. Both of these applications essen- 
tially require knowledge of the total Galactocentric solar 
velocity, V^©, to transform velocities into the Galacto- 
centric (dark- matter-halo) rest frame 25 . Although this 
motion is often decomposed as V c plus the Sun's mo- 
tion with respect to V c , our results in this paper show 
that this is unnecessarily dangerous, as it requires the 
strong assumption that the LSR is on a circular orbit, 
such that V^© = V c + V^q^lsr, where the final term 
is the solar velocity with respect to the LSR. Our re- 
sults in § 4 strongly rule out that this assumption holds 
for the currently accepted solar motion of ~ 12 km s" 1 . 

25 The one exception is that in the Standard Halo Model, the 
velocity dispersion of the isothermal halo is given by V c /y/2. How- 
ever, the velocity dispersion of the halo is better determined di- 
rectly from observations than by using this assumption. 
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However, this approach is entirely unnecessary for cor- 
recting velocities to the Galactocentric rest frame, as the 
proper motion of Sgr A* combined with the "clean" esti- 
mate of Rq from Galactic-center dynamics shows that 
V^o ~ 242 km s" 1 , which agrees with our measured 

value of V^q. Adopting a standard V^q = 242 km s _1 
for correcting velocities for the Sun's motion, and decou- 
pling this correction from the question of the true value 
of V c , would therefore be helpful. 

5.4. Implications for the Mass of the Milky Way 

We can use our measurement of the inner rotation 
curve of the Milky Way to estimate the Milky Way's total 
dark-halo mass. We combine the measurements in this 
paper with those of the Milky Way's outer rotation curve 
(R > 20 kpc) of Xue et al. (2008) (specifically, taking 
the measurements using simulation II in Table 3 of that 
paper). Based on FIG. 6 and TABLE 2, we add our mea- 
surement as a flat rotation curve with V c = 218 km s - , 
evaluated at three radii with uncertainties of 3 km s _1 , 
6kms _1 , and 10 km s -1 at R = 4 kpc, R = 8 kpc, 
and R = 12 kpc, respectively (removing the measure- 
ments at R = 7.5 kpc and R = 12.5 kpc from Xue et al. 
2008). We also use our recent measurements of the local 
density of dark matter, p DM = 0.008 ± 0.003 M pc~ 3 
(Bovy & Tremaine 2012), and of the disk scale length, 
h R = 3.25 ± 0.25 kpc (Bovy et al. 2012b), revising the 
measurement from the latter paper downward to correct 
for the influence of the dark halo on that measurement. 

We model the Milky Way's potential as the combi- 
nation of a bulge component with a Hcrnquist profile 
with a scale length of 600 pc, a disk component with 
a scale height of 300 pc (Bovy et al. 2012a) and a scale 
length that is a free parameter, and a Navarro-Frenk- 
White halo (NFW; Navarro et al. 1997), with a scale ra- 
dius that is a free parameter (i.e., without constraint on 
the concentration). The relative contributions of these 
three components are free parameters. We then fit the 
data given in the previous paragraph, and we find that 
Mhaio = x 10 11 Mq. The total mass of the Milky 
Way is -8.5 x lO n M . 

This "low" estimate for the mass of the Milky Way 
likely makes it less massive than M31, which has an es- 
timated dark-halo mass of 14 x 10 11 Mq (Watkins et al. 
2010). This measurement leads to a combined mass for 
the Local Group of —2.4 x 10 12 M Q , which is consistent 
with the Local-Group-timing-argument within the large 
cosmic scatter (van der Marel et al. 2012). Our measure- 
ment of V c = 2 18 ±6 km s~ combined with its estimated 
/-band magnitude of —22.3 mag makes the Milky Way 
undcrluminous with respect to the Tully-Fisher relation 
of external spiral galaxies by about la (see Flynn et al. 
2006). 

6. CONCLUSION 

In this paper, we have measured the Milky Way's ro- 
tation curve over the range 4 kpc < R < 14 kpc from 
the new APOGEE data set of kinematically-warm stel- 
lar tracers at large distances from the Sun. Our mea- 
surement is not "clean" , in the sense of being a geomet- 
ric measurement such as that provided by tangent-point 
observations of HI emission or by a measurement of the 



Sun's motion relative to an object assumed to be at rest 
with respect to the Galactic center. Because we use a 
warm stellar population, we must correct for the offset 
between the average rotational velocity of this population 
and the circular velocity — the asymmetric drift — using 
a dynamical Jeans model. However, the fact that our 
measurement uses a dynamical effect, rather than being 
purely geometric, has the advantage that we unambigu- 
ously measure V c (R) in the sense of the radial force at R, 
and that we can measure the solid-body-rotation contri- 
bution to the rotation curve, in contrast to HI measure- 
ments. 

Our main results are discussed in § 4. We find that the 
Milky Way's rotation curve is approximately flat over 
4 kpc < R < 14 kpc, with V C (R ) = 218 ± 6 km s _1 . 
TABLE 2 summarizes our results, and provides some 
alternative representations, such as the Oort constants, 
the local rotational frequency, and the contribution of 
the non-flatness of the rotation curve to a cylindrical- 
Poisson-equation determination of the local surface-mass 
and dark-matter densities. We simultaneously measure 
the Sun's velocity in the Galactocentric rest frame — 
these are independent free parameters in our model — and 
find that Vr,© = — 10±1 km s _1 and the angular motion 
of the Galactic center fi = 6.3+g'7inas yr _1 . These val- 
ues are consistent with previous studies. Our measure- 
ment of = 242^3° km s _1 then leads to a solar offset 
from V c that is larger than the locally-measured value 
(~ 12 km s" 1 ; Schonrich et al. 2010) by 14 ± 3 km s _1 . 
This result may indicate that the challenging local mea- 
surement of the solar motion is incorrect or that the so- 
lar velocity is influenced by non-axisymmetric streaming 
motions that put the LSR on a non-circular orbit. 

Looking forward, we expect to accurately measure 
the atmospheric parameters and abundances other than 
[Fc/H] from the high-resolution APOGEE spectra in the 
near future. This will improve the measurements de- 
scribed in this paper by (a) providing better dwarf/giant 
separation, (b) allowing for more direct and more precise 
distances to be derived, and (c) letting the analysis be 
conducted on different chemically-defined populations of 
stars. In particular, this will allow for a more detailed 
investigation into the signatures of non-axisymmetric dy- 
namics in our data. 

The data upon which the measurements presented in 
this paper are based will be released as part of SDSS-III's 
Data Release 10 in the summer of 2013. 



It is a pleasure to thank Scott Tremaine for many illu- 
minating discussions about this work. We also thank 
Jennifer Johnson, Mark Rcid, Hans- Walter Rix, and 
Greg Stinson for helpful comments. J.B. was supported 
by NASA through Hubble Fellowship grant HST-HF- 
51285.01 from the Space Telescope Science Institute, 
which is operated by the Association of Universities 
for Research in Astronomy, Incorporated, under NASA 
contract NAS5-26555. J.B. was partially supported by 
SFB 881 funded by the German Research Foundation 
DFG and is grateful to the Max-Planck Institut fur As- 
tronomic for its hospitality during part of the period dur- 
ing which this research was performed. T.C.B. acknowl- 
edges partial support by grants PHY 02-16783 and PHY 



19 



08-22648: Physics Frontiers Center/Joint Institute for 
Nuclear Astrophysics (JINA), awarded by the U.S. Na- 
tional Science Foundation. This research made use of the 
emcee MCMC sampler (Foreman-Mackey et al. 2012). 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 

This publication makes use of data products from the 
Wide-field Infrared Survey Explorer, which is a joint 
project of the University of California, Los Angeles, 
and the Jet Propulsion Laboratory /California Institute 
of Technology, funded by the National Aeronautics and 
Space Administration. 

Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, and the U.S. Department 



of Energy Office of Science. The SDSS-III web site is 
http: / /www. sdss3.org/. 

SDSS-III is managed by the Astrophysical Research 
Consortium for the Participating Institutions of the 
SDSS-III Collaboration including the University of Ari- 
zona, the Brazilian Participation Group, Brookhaven Na- 
tional Laboratory, University of Cambridge, Carnegie 
Mellon University, University of Florida, the French 
Participation Group, the German Participation Group, 
Harvard University, the Instituto de Astrofisica de Ca- 
narias, the Michigan State/Notre Dame/JINA Participa- 
tion Group, Johns Hopkins University, Lawrence Berke- 
ley National Laboratory, Max Planck Institute for As- 
trophysics, Max Planck Institute for Extraterrestrial 
Physics, New Mexico State University, New York Uni- 
versity, Ohio State University, Pennsylvania State Uni- 
versity, University of Portsmouth, Princeton University, 
the Spanish Participation Group, University of Tokyo, 
University of Utah, Vanderbilt University, University of 
Virginia, University of Washington, and Yale University. 



APPENDIX 

PHOTOMETRIC DISTANCE PDFS FOR GIANT STARS 

In this appendix, we describe the calculation of the photometric-distance distributions p(d\l,b,(J — 
K s )o, Hq, [Fe/H], DF, iso) for the stars in our sample. The photometric-distance PDF for each star is the combi- 
nation of the probability of its observed photometry, given a model for the color-absolute-magnitudc distribution of 
giant stars, and a prior on the distance. Thus, we write 

p(d\l, b, ( J - K s ) , H , [Fc/H], DF, iso) <x p((J - K s ) , H \d, [Fc/H], iso) p(d, 1, 6|DF) , a x 

[ Al ) 

= iso[ Fc/H ](-ffo - li{d), (J - K a )o) v*(R,z\DF) d 2 cosb, 

where /i is the distance modulus, and d 2 cosb is the Jacobian for the coordinate transformation (R,<j),z) — > (d,l,b), 
because we have anticipated that the distance prior is written in Galactocentric cylindrical coordinates (R,z,(j)). This 
prior is an exponential distribution in R, v se (R, z|DF) oc exp (—R/hjt); our sample is sufficiently close to the plane that 
vertical density gradients are unimportant (see FIG. 2). We have implicitly assumed that the prior does not depend 
on [Fc/H]. Our fiducial model has Iir = 3 kpc, as is appropriate for the metal- rich disk stars that make up our sample 
(Bovy et al. 2012c). 

In Equation (Al), we have also anticipated that we obtain the probability of the observed photometry of the star, 
given its distance from a model for the isochrones of giants, combined with a model for the distribution along the 
isochrone from an initial-mass-function model; this probability is labeled as 'iso'. We only use the color (J — K s )q and 
the magnitude H , as little is gained by adding J or K Sj0 separately, and we compute the density in the color-absolute- 
magnitude plane for a given age and metallicity (which constitutes a single isochrone) by counting the number of stars 
generated from an initial mass function in small boxes in (J — K s )q - Hq. For stars in our sample, we find the nearest 
Z/Zq (assuming Z & = 0.019) on a grid with a 0.0005 dex spacing, and we average Padova isochrones for metallicities 
Z/Zq, Z/Zq - 0.0005, and Z/Z Q , +0.0005, except at the edges of the grid at Z/Z Q = 0.0005 and Z/Z Q = 0.03. We 
average over age by assuming a constant star-formation rate up to 12 Gyr (except for the fit in § 4.2 with multiple 
populations, where we use an exponentially-declining star- formation rate). We use a lognormal Chabrier (2001) model 
for the IMF (but the resulting distance PDFS are essentially the same when using a Kroupa 2003 IMF). The resulting 
density for a star with solar metallicity is shown in FIG. 11. For a given distance and [Fe/H], the probability of a 
star's magnitude given its color is evaluated by computing its absolute magnitude Mjj and evaluating a density like 
the one shown in FIG. 11. 

The distribution isope/H] (-^o — /- l (d),(J — K s )q) is relatively insensitive to changes in [Fe/H] at high metallicity 
([Fe/H] > —0.4), where 90% of our sample lies. Thus, for the majority of the stars in our sample, the photometric 
distance PDF does not strongly depend on the measured [Fe/H] and resembles FIG. 11. 

Also shown in FIG. 11 is the direction that extinction moves objects in color and magnitude. In contrast to the 
case for the main sequence, this direction is at a large angle with the locus of the giant branch, such that under- or 
over-estimated extinction values can significantly change the photometric-distance PDF. The length of the extinction 
arrow in this figure is the median extinction in (J — K s ) for the sample (0.45 mag). 

ANALYSIS TESTS ON MOCK DATA SAMPLES 

In this section, wc perform extensive mock data analyses to test the methodology and approximations described 
in § 3, and to determine at what level the data are sensitive to changes in the Galactic parameters. We cre- 
ate mock data sets by re-sampling the line-of-sight velocity for each data point from its PDF p(V\ os \l, b, (J — 
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Fig. 11. — Density in the color— magnitude plane, for the giant branch of stars with solar metallicity. Each color bin is normalized 
separately to show the absolute-magnitude distribution at each color. This density is calculated using Padova isochroncs (Bertelli et al. 
1994; Bonatto et al. 2004; Marigo et al. 2008; Girardi et al. 2010), assuming a lognormal Chabrier (2001) IMF and a constant star-formation 
rate. This figure is an example of the density that we use as the photometric-distance distribution for each star. The length of the extinction 
arrow is the median extinction for the sample (0.45 in J — K B ). 



TABLE 3 

Galactic Parameters for a Flat-Rotation-Curve Fit to Mock Data Sets 



Parameter Data / Mock input Mock 1 Mock 2 Mock 3 Mock 4 Mock 5 



V c (Ro) [kms- 1 ] 

Ro [kpc] 

V R , e [km s- 1 ] 

Qq [km s kpc ] 

vr(Ro) [km s _1 ] 

Ro/ha 

y2 — _2 /_2 


218±6 

-10.5t ; 6 8 
30.0+3 3 
31.4t°i 
03 +0 1 

o 7n+ - 30 

u -' u -0.01 


225l 10 

-10.3±g-2 
30.7l°j 

32.0+22 
o in+° :o1 

n OO+0.21 


222±5 

-9-3±g:S 

30.6+3 2 

28.8l ;° 

-0.07±0.06 

60+ 012 
U - DU -0.03 


221 + 12 

7.91H 
30.3+°* 

3i.5+i ; ; 

o n7+"" 4 

u - u '-0.09 

n s9+° 13 


225+^ 

8.01°;? 

-io.8±8:S 

30.7 + ^| 

34.3+3 ,: i 

25+ ' 01 
UiZO -0.19 

35+° ' 19 


217li 2 

7 q+0.9 

'■°-o.i 
-9 4+°.'? 

0.4 

30.31°;* 
29.3lJ'| 
_n n4+ OS 

u - u *-0.04 

o 5n+ 09 

U.OU_ Q3 


Ax 2 /dof 




-0.19 


-0.18 


-0.20 


-0.17 


-0.22 



Note. — The final line of this table has the difference in x 2 per degree-of-freedom between the best-fit to the 
mock data set and the best-fit to the real data. Parameters are as in TABLE 2. 



K s )o, Hq, [Fc/H], V c (R), Ro, Vr ; 0, V^©, DF, iso) (Equation (1)) in the best-fit model with a flat rotation curve for the 
real data (TABLE 2). However, rather than using the simple Gaussian-with-asymmctric-drift-offsct model, we sample 
from a Dehnen distribution function (Equation (6)), with a radial scale length of Iir = 3 kpc, and a radial- velocity- 
dispersion scale length taken from the best-fit model. Because the Dehnen distribution function has X 2 built-in, we 
do not use the best-fit X 2 . For a Dehnen DF, the marginalization of the planar velocity distribution over the com- 
ponent of the velocity that is tangential to the line-of-sight cannot be done analytically, and we numerically integrate 
over this component. We do not add any observational error to the mock line-of-sight velocities, as the observational 
uncertainties are vanishingly small (see § 2). In this manner we produce five mock data sets that are exactly like the 
real data, except for the line-of-sight velocities. 

We apply the exact same fitting procedure to the mock data as is used for the real data. We fit the five mock data 
sets with a flat-rotation curve model, and obtain the best-fits and uncertainties given in TABLE 3 (for comparison, 
this table also includes the best-fit model for the real data that was used to generate the mock data). We see that 
the true values for the parameters of interest are recovered well by the methodology of § 3. We also see that the 
uncertainty ranges on the Galactic parameters are similar, both in size and asymmetry, to those derived for the real 
data. Therefore, the best-fit parameters do not appear to be biased by the approximate methodology used in this 
paper, and the uncertainties arc as should be expected from these mock data tests. 

We also fit models with a non-flat rotation curve to the mock data sets. We only do this for three of the mock data 
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TABLE 4 

Galactic Parameters for non-Flat-Rotation-Curve Fits to Mock Data Sets 



Rotation curve 



Parameter 



Data / Mock input Mock 1 



Mock 2 



Mock 3 



power-law 



V c (Ro) [km s- 1 ] 


Ro [kpc] 



R,r._, 



[km s 



£1q [km s 1 kpc 1 ] 
a R (Ro) [km s- 1 ] 

X 2 



218±6 



.V 



-1.2 
-0.1 

-10.5+ 05 



-0.8 
1+0.3 

31.4+°'i 



30. i 



-3.2 
,+0.01 
-0.27 

70+ - 30 
u -' u -0.01 



0.03 i 



224 



+28 
2 



07+ - 06 



1+1.2 
-0.1 
7+0.7 
-0.4 

30.8+Jj 

j+0.2 
'-2.2 
7+0.01 
-0.14 
,,+0.16 
-0.01 



8.0: 
-9.7" 1 



31. 

0.071 

0.44 4 



,+8 

-21 

02+ - 02 
u - uz -0.11 
1+0.7 
-0.1 
5+0.5 
-0.7 
4+0.2 
-3.9 

29.4±1.0 

-0 n5+ a05 
U - Ua -0.08 

59 +0 10 
u ' o9 -0.04 



222" 1 



8.0Z 
-9.3j 
30.4"* 



H4 
-27 
1-0.01 
-0.14 
1+0.7 
-0.1 
_Q 1+0-7 
a ' x -0.6 

33.3l°* 
7 +o!o4 

-0.11 

36+ 09 
u " 3D -0.02 



224] 
0.08j 
8.0^ 



0.17"* 



linear 



V c (Ro) [km s- 1 ] 

dVc/dR (Rq) [km s" 1 kpc- 1 ] 

Ro [kpc] 

V RtQ [km s- 1 ] 

Hq [km s _1 kpc ] 

a R (Ro) [km s- 1 ] 

Ro/h a 

Y 2 — r, 2 Irfi 



21S±6 



1+1-2 
1 -0.1 



-10. 5 -1 



:;+0.5 
-0.8 

30.0+ 03 



-3.3 
4+0.1 

0.03+°' §1 



31. 4 -1 



0.70 



0.27 
+0.30 
0.01 



223±12 
2.0+.°;? 

8-oto:? 

-9.7±0.7 

R +0.3 
'-4.3 

31.3l°-| 

4+0.03 
-0.09 
, + 0.13 
-0.03 



30.1 



0.04: 

0.46 4 



222" 1 



5+20 
1+1 

8.0+1.? 

-9-i±8:l 

31+ 1 



30.0 
-0.03 
0.59 



3 

+0.2 
2.1 
+ 0.01 
0.13 
+ 0.14 
0.01 



2.8±gl 

8 0+ o:s 
°- u -o.i 

-8.9±0.6 

30.6i^| 

32.8l5 : f 

12+° : ° 2 
u ' lz -0.10 



0.40 



+0.13 



cubic 



V c (Ro) [km s- 1 ] 
dV c /dR(Ro) [km s- 1 kpc -1 ; 
d 2 V c /dR 2 (R ) [km s- 1 kpc" 



218+6 



d 3 V c /dR 3 (Ro) [kr 
Ro [kpc] 



R..-. 



[km s 



£1q [km s 1 kpc ] 
ur(Ro) [km s- 1 ] 
Ro/h a 



kpc 



.r 



1+1.2 
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-0.8 
1+0.3 

31.4+° :i 



30. 0^ 



-3.2 
,+0.01 
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u -' u -0.01 



0.03- 1 



,+21 
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002+ 040 
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Parameters in this table are as in TABLE 2. 




Ra < 9 kpc 
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Fig. 12. — Same as FIG. 6, but for mock data set 2. To completely mimic the procedure applied to the real data we have imposed a 
Ro < 9 kpc prior for the cubic-polynomial fit to the rotation curve, even though this only excludes tens of the 10,000 points that were 
sampled from the PDF (see TABLE 4). 



sets; the results for fitting a power-law rotation curve, or a linear or cubic polynomial curve, are given in TABLE 4. 
We see that, in these cases, the Galactic parameters are recovered within the uncertainties, with perhaps a slight bias 
of a few km s _1 in the best-fit parameters, although the (asymmetric) PDFs for the best-fit parameters do contain the 
"true" values within the 68% confidence intervals given. Again, the size of the uncertainties on the best-fit parameters 
for the mock data are similar to those found for the real data in TABLE 2. 

In FIG. 12, we show for one of the mock data sets the same representation of the constraints on the shape of the 
rotation curve as is given for the real data in FIG. 6. Comparing these two figures, we see that the constraints are 
similar, with the real data giving a slightly narrower range around the best-fit, which is approximately fiat. It is clear 
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that around R = 14 kpc, the constraints on the rotation curve from the current APOGEE data become weak, so we 
limit our discussion to 4 kpc < R < 14 kpc in this paper. 
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